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A method to calculate the mean squared matrix element of weak interaction between compound 
states is developed. The result is expressed in terms of matrix elements of the nucleon-nucleon strong 
and weak interactions times the Fermi distribution functions at finite temperature. Numerical cal¬ 
culations for ^3®Th are in excellent agreement with recent measurements of parity nonconservation 
effects in neutron capture. In fact, our calculations prove that the factor of dynamical enhance¬ 
ment (ratio of compound-nucleus effect to single-particle one) really exceeds 100, thus making it 
unnesessary to assume a value of the weak constant bigger than standard one {g ~ 10® e ~ 1-^4). 
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The parity nonconcerving (PNC) nucleon interaction 
in nuclei now attracts attention of both experimental¬ 
ists and theorists, especially, in connection with the re¬ 
cent experiments on slow neutron scattering through the 
compound nuclear states, where the measured PNC effect 
(dependence of cross section on neutron helicity) proved 
to be of order several per cent i, i (cf. with the PNC 
effects in p-p or p-a scattering where effect is of order of 
3 • 10“^). Moreover, correlation of sign in the effect on 
close neutron resonances has been observed |^. 

In the current literature |3|-||] , two different approaches 
to explain a great value of this effect coexist, based on 
different assumptions and contradicting each other. The 
first one, the statistical model of dynamical enhancement, 
was considered in works i- 0 .(In fact, a large value of 
the effect was predicted in the I980-I98I papers of Ref. 
0.) Within this approach, the essential enhancement of 
the parity violating amplitude in neutron capture arises 
from the mixing of closely lying (within interval of sev¬ 
eral eV) compound nuclear levels of opposite parity, and 
statistical enhancement (~ ^fN) of the weak matrix el¬ 
ement between compound wave functions composed of 
N ^ 3 X 10® many-particle configurations. The estimate 
of the magnitude of the effect was given in §, i, i, 
based on the standard Hamiltonian of the weak interac¬ 
tion of nucleons in a nucleus 

W =-^-{(crp),p}, e= 1 . 0 - 10 -® 5 , (1) 

2v 2m 

where G = is the Fermi constant, m is the 

nucleon mass, p and a are the neutron momentum and 
its doubled spin correspondingly, while p is the nuclear 
density; the nucleon dimensionless constant Pp^n (see e.g. 
Ref. H) is of order unity (now the notation e is widely 
used Q). 

The second approach, the so called “valence mecha¬ 
nism” ||] (see also |^), is based on the assumption that 
the weak amplitude admixing the s-wave to initial p-wave 
is dominated by direct matrix elements of PNC-potential 
(1) between these states, so the effect is assumed to be 


of single-particle nature. The valence mechanism gives 
correlations of the sign in the effect on the different reso¬ 
nances but to explain the observed magnitude of the ef¬ 
fect in this approach, one has to use the neutron constant 
in Eq.(l) being at least two orders of magnitude larger 
than is predicted by theory (see Refs. |^, [^). Thus, 
the two mechanisms require a neutron weak interaction 
constant which differs by two orders of magnitude in the 
range 1-300. It becomes even more important in connec¬ 
tion with the so-called “Tsinoev puzzle” where the 
observed PNC effect is 10® times bigger than the theo¬ 
retical estimate. 

Staying within the framework of the first approach (dy¬ 
namical enhancement) we present here a method to cal¬ 
culate the mean squared weak matrix element between 
s and p compound resonances, in the statistical model 
with account for nuclear structure and realistic residual 
nucleon interaction. We should note that the matrix ele¬ 
ment between compound states was considered earlier in 
Refs. However, these works use some hypotheses 

which are not easy to justify. In the work Q a pro¬ 
portionality relation between the matrix elements of the 
weak and residual interaction was used. In the work jj] 
it was supposed that the matrix element between com¬ 
pound states is given by the same formula as the matrix 
element of the nucleon excitation from the ground state 
with only some minor modifications (occupation num¬ 
bers and the optical potential depend on the tempera¬ 
ture of the compound nucleus). In our approach we have 
not used these hypotheses, and our result looks different 
(e.g., it is not proportional to the square of the frequency 
of the time-dependent field or T® (w = 0 for the weak 
interaction, and T is a temperature)). 

Note, that we do not discuss sign correlations in the 
present paper. 

Calculation of mean squared weak matrix element is 
based on the equivalence theorem of canonical and mi- 
crocanonical ensembles for a system with a large number 
of degrees of freedom. Let us remember that the wave 
function of any compound state with angular momentum 
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j and parity tt may be expressed as follows 
\f)='^Ca\a>, |a >= (a+5c+cie+...)j.r|0 >, (2) 

Oi 


through their components \a > being many-particle ex¬ 
citations over the shell-model ground state |0 > (we will 
denote them by simple Dirac brackets saving the nota¬ 
tion I) for compound wave functions; here and in what 
follows notation means the coupling of nucleon cre¬ 

ators o'*" and destructors a to total angular momentum 
j and parity tt) . Amongst them, it is reasonable to sepa¬ 
rate explicitly the contribution of the “principal compo¬ 
nents”, dominating normalisation of the compound 

state, Eq.(2). The energies of these components must be 
within the interval [E — where E is the 

energy of the compound state and Pspr is the spreading 
width of the component (typically, P^jp^ ^ 2MeV, Refs. 

These components (which contain several excited 
nucleons) can be composed by excitations of protons and 
neutrons only inside incomplete (valence) shells. Mean 
squared values of the coefficients C{Ea can be described 
by the formulae (see e.g.. Ref. [p^) 


_ 1 _ ttT^ 

CHE^) = =A(T,pr,E-E^), N = 


( 3 ) 


A(rspr,E - Ea) = 


/4 

spr / 


iE-Ea)^+T%jA^ 


where Ea is the energy of an arbitrary many-particle con¬ 
figuration), d is the averaged energy distance between 
the resonances, and N is the number of principal compo¬ 
nents. The Breit-Wigner-type factor A, describing cut¬ 
ting off of weights before states distanced in energy, may 
be treated as a “spread” J-function normalized as to be of 
order unity for \E — Ea\ < Tspr/'^ and with conventional 
limit A{Tspr,E — Ea) ^^^’”' S{E — Ea) for Tsp^. ^ 0. 
Thus, the coefficients before the “principal” components 
Ca in (2) are governed by the microcanonical ensemble 
rule |12|,[7|. The very important point should be born in 
mind, that there are no single-particle states of opposite 
parity having the same angular momenta within the va¬ 
lence shells. Since the weak interaction (Eq.(l)) mixes 
only states of such type, it follows from the definition 
of “principal” components that the weak matrix element 
between two compound state of close energy is dominated 
by the weak transitions between “small” components of 
one resonance and “principal” ones of the second reso¬ 
nance, and vice versa (this was firstly mentioned by the 
authors of Ref. Q, see also 0). Any transfer of one par¬ 
ticle from the valence shell to another one gives a rise 
in excitation energy as large as Egp ~ 8MeV (what is 
much more than typical matrics elements of the resid¬ 
ual interaction V) leading out a configuration from the 
microcanonical ensemble of “principal” components ac¬ 
cording to Eq.(3). Therefore, one can easily generate 
the appropriate set of “small” configurations by means 
of first-order perturbation theory in the residual strong 


interaction V. Thus, matrix elements of the weak inter¬ 
action between compound states look like 

E Ea 

a 

^ ((s|lU|/3 X/3|Ub)) 

^ E-Eg ’ ^ ^ 

where \a > and |/3 > are small components, and com¬ 
pound states |s)). Ip)) contain only principal components. 
We stress that we do not need any “exotic” parts of 
the residual strong interaction here. Since E — Ea ~ 
8MeV ^ V (see above), only the dominating and well- 
known parts |13|-|l5| of the two-nucleon interaction 

V = ^ X! a'^bVab,cdC'^d, (5) 

ab,cd 


which will be specifyied below, are important in 
Eqs.(4),(5). Here, our consideration is general and not 
even confined to nuclear system. 

The weak interaction (Eq.(l)) is a single-particle op¬ 
erator. This fact allows one to the include weak inter¬ 
action into the mean nuclear field and transfer the per¬ 
turbation theory expansion in the single-particle orbitals: 
i’a = '>Pa + La where Ca and ca are the 

energies of the orbitals ipa and tj^A (differing by their par¬ 
ities), the large Latin indices label the corresponding off- 
valence-shell states. Thus, we can express result in terms 
of the residual interaction renormalised by the weak in¬ 
teraction {Vab,cd = V{db,cd)): 




abed 




VAb,cd 
Ca - CA 


-WaA 


E 


^aB ,cd 

Sb — e_B 


-WBb 


E Vab,Cd 

- - -w,c 


Cc - ec 


E Vab,cD 

- WDd, 

D 


( 6 ) 


here WaA =< V'A|lU|'!/'a >. With these notations, to 
the first order in V, Eq.(4) can be read as follows: 
{s\W\p) = ((s|U|p)). The advantage of using the ef¬ 
fective two-particle PNC-interaction V is that the ma¬ 
trix elements between compound states are expressed 
through the matrix elements Vab,cd between valence shell 
single-particle states. Thus we avoid the necessity of ex¬ 
plicitly considering the “small” components of the com¬ 
pound states which we believe cannot be described by 
the same spreading widths as the principal components 
(see Eq.(3)). 

Consider now the mean squared value of this matrix 
element: 


{p\W\s)[s\W\p) = ((p|U|s))((s|U|p)) = 

= Y,CaCg{{p\V\a><l3\V\p)) = 

oi(3 

= E E - Eo.)mV\a X a\V\p)). (7) 
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Here, we have expanded the compound state |s) in terms 
of the components (2) and made use of the statistical 
independence of the coefficients Ca (see Eqs.(2,3),Ref. 

§, §, i): 

( 8 ) 

In the second quantization representation, summation 
over a in (7) is equivalent to summation over different 
components of operator the V in Eq.(5), i.e. the prob¬ 
lem is reduced to the calculation of {{p\VV\p)). Then, 
to calculate the averaging over p-resonance “principal” 
components ((p|...|p)) in let us use, instead of the 
present microcanonical ensemble, the equivalent canon¬ 
ical one (which may be choosen for a system with a 
large number degrees of freedom by introducing the ef¬ 
fective nuclear temperature T and chemical potentials 
A„,Ap). In such a way, the expectation value in (8) is 
reduced to a canonical ensemble average with the stan¬ 
dard contractor rules {{p\a^b\p)) = Sab^a ^ for zzj being 
the finite temperature Eermi occupation probabilities, 
= {exp[{ea — A)/T] -|- 1}“^. The canonical ensem¬ 
ble parameters T, (r means isospin projection) are to 
be determined from conventional “consistency” equations 
E = Y,a Z = J2p and N = J2n for the exci¬ 
tation energy E (being equal the to neutron separation 
energy, Bn), nuclear charge Z, and neutron number N 
correspondingly. After contractor evaluations one simply 
obtains from (7): 

(9) 

y 'Tri spr 

^ (1 Vab.cd - Vad,cb P 

^ abed 

1 ^ 

A(rspy., Ca Cfo -f Cc Cd) / . (10) 

The argument of the function A here is the change of the 
energy: E - E^ = Ca - Cb + Cc - Cd, and V is given by 
Eq.(6). In fact, it is an approximate energy conservation 
law with the accuracy up to width of states. 

The numerical calculations for ^^^Th have been per¬ 
formed with the use of single-particle basis of states 
obtained by numerical calculations of the eigenvalue 
problem for the Woods-Saxon potential with spin- 
orbital interaction in the form U{r) = —Uof{r) + 
Uis{<7\){h/+ Uc with /(r) = (1 -f exp{{r - 
R)/a))~^, where 1 is the orbital angular momentum, 
Uc means Coulomb correction for protons, Uc = 
32'e^/(2i?)(l — r^/(3i?^)), r < R and Uc = Ze^ jr, r > R, 
for R, a, and r being the nuclear radius, diffusity param¬ 
eter and radial variable correspondingly. The parame¬ 
ter values were used in accordance with Bohr-Mottelson 
formulae (see Ref. [Q) for the case of they are 

close to those established for heavy nuclei like lead (Ref. 
0 ) to reproduce single-particle properties. As for the 


residual interaction, we have employed the most widely 
used Landau-Migdal particle-hole interaction of contact 
type with spin- and isospin-exchange terms which rises 
to Landau Fermi liquid theory (Ref. Q); for the case of 
a nucleus it was established in the Theory of Finite Fermi 
Systems 00 by summation of all graphs irreducible in 
the particle-hole direction. In that case, the explicit form 
of the matrix V in (5) is given by the second quantized 
version of the interaction 

U(r r') = CS{r — r')[f + J'tt' -\- gerer' + g'Tr'aa'], (11) 

where C = 300 MeV fmP is the universal Migdal con¬ 
stant 0, 0, tT(cr') and t { t ') mean particle(hole) spin 
and isospin Pauli matrices, respectively, and the strenghs 
f,f',g,g' are in fact functions of r via density depen¬ 
dence: f = fin- {fex - fm){p(r) - p(0))/p(0) (the same 
for g,g') with ansatz p{r) = p(0)/(r) (see above). (Quan¬ 
tities subscripted by “in” and “ex” characterize inter¬ 
action strenghs in the depth of the nucleus and on its 
surface, respectively). This interaction, with its param¬ 
eter values listed below, has been successfully used by 
many authors (see Refs. [0) to quantitatively describe a 
great amount of various properties of heavy nuclei. The 
value of temperature T = 0.6MeV was used in accor¬ 
dance to the consistency condition for excitation energy 
(see above). 

We present here the results for the conventional choice 
of values for parameters in the Landau-Migdal interac¬ 
tion, Eq.(lO), which has been widely used for heavy 
nuclei (see 0, 0 and references therein), namely 
fex = -1.95, /„ = -0.075 = 0.05 fl^ = 0.675, 

gin = gex = 0.575, and g[^ = g'ca; = 0.725. Note that 
the exchange matrix elements of V enter in Eq.(9) but, 
generally, the values of the parameters /, /', g, and g' 
are chosen in such a way that the exchange terms are al¬ 
ready included (it can always be done by use of the Firtz 
transformation). The variant referred to below as I cor¬ 
responds to this standard procedure; we present also the 
results for the case when the excange terms in (9) are in¬ 
cluded explicitly (referred to as II). The results for \/w’^ 

(8) may be expressed explicitly in terms of the proton 
and neutron weak constants gp and gn in the form: 

^/W^=^^/W^{Eppgp)^ + {E 

nnQn )2 -I- EpngpgnimeV); 

( 12 ) 

the numerical calculations of the constants E gives the 
following results: 

I : = 2.08meV; II : VW = 3.57meU. 

The experimental value is ( (H). 

The values of \/w^ for the cases I and II were obtained 
for the conventional choice of weak constants, gp = A and 
= I (Ref. 01). (In the notation widely adopted in the 
current literature [H, for the neutron case it corresponds 
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to the value e = 10“®5„ = 10“®). Now, we can com¬ 
pare this result with single-particle (valence) estimation 
Wyai ■ In the valence mechanism, only single-particle com¬ 
ponents contribute (in ^®®Th, 4s and Ap neutron states). 
Therefore, the valence model result is 

w^ai - = < 4s|lU|4p >~ ^5„0.740(eU) = 1.72 • 10-®meU. 

(13) 

Thus, statistical contribution is 10® times bigger (due 
to the extra factor VW, compare Eq.(9,ll) with (12)). 
This factor reflects the incoherent contribution of all 
N components. Calculations fulfilled for other sets of 
parameters display no strong sensitivity of the numeri¬ 
cal results to variations of both single-particle basis and 
residual interaction. Let us point out that in our ap¬ 
proach, the only essential assumption made is that of 
the statistical properties of the distribution for coeffi¬ 
cients Ca in Eqs.(2),(8). As uncertainity in definition 

of N is concerned, two estimates for it, namely 
N ~ ^ ~ 17eU in 

spreading T^jp^ — 2MeV) and that from the widths of 
the compound resonances, N ~ give 

approximately equal values N ~ 2.3 • 10“® with ac¬ 
curacy up to a factor of 2 (here T® is the width of the 
s or p single-particle resonance (see Ref. 0 ) , and r„ is 
the width of the s or p compound resonance correspond¬ 
ingly). 

The results of this work can be summarised as follows. 

A consistent method to calculate the mean squared weak 
matrix element between two compound states of oppo¬ 
site parity is proposed, based on a statistical model with 
account for nuclear structure. The results prove that the 
dynamical enhancement ~ ViV does really exist. The ob¬ 
served large value of this quantity for ^®®T/i is explained 
in terms of the model with the conventional neutron weak 
constant. Moreover, even for the zeroth value of the neu¬ 
tron constant the effect is still reproduced without sub¬ 
stantial cut-off because of the proton PNC transitions in 
the nuclear compound state. (In fact, even for = I, 
proton contribution is few times larger due to the bigger 
value of constant gp). Further experiments in this region 
would be of great importance. 

We are grateful to V.F. Dmitriev and V.B. Telitsin for 
kindly providing us with a computer program to solve 
eigenvalue problem in Saxon-Woods potential. 
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